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Parquet equations, describing the competition between superconducting and 
density-wave instabilities, are solved for a three-dimensional isotropic metal in a high 
magnetic field when only the lowest Landau level is filled. In the case of a repulsive 
interaction between electrons, a phase transition to the density- wave state is found at 
finite temperature. In the opposite case of attractive interaction, no phase transition 
is found. With decreasing temperature T, the effective vertex of interaction between 
electrons renormalizes toward a one-dimensional limit in a self-similar way with the 
characteristic length (transverse to the magnetic field) decreasing as \nr'^^^{ijJc/T) 
{u)c is a cutoff). Correlation functions have new forms, previously unknown for con- 
ventional one-dimensional or three-dimensional Fermi-liquids. 
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The behavior of an isotropic three-dimensional (3D) metal in a high magnetic field has 
attracted attention of physicists for a long time. In this system, the energy of an electron 
depends only on the momentum along the magnetic field. Thus, the system exhibits effects 
characteristic of one-dimensional (ID) metals, while intrinsically it is 3D. As an immediate 
consequence of this fact, it was suggested that the system should be unstable with respect 
to charge- or spin-density- wave (DW) formation [|l|. Another suggestion was that the sys- 
tem could remain superconducting (SC) in an arbitrarily high magnetic field, since the ID 
dispersion law still allows for a SC instability 0. It was pointed out already in Ref. that 
both instabilities must be taken into account simultaneously in the so-called parquet approx- 
imation. The parquet equations were written correctly in Ref. 0, where it was found that 
the DW solution is indeed an asymptotic solution of the equations. However, the equations 
were not solved numerically, thus the question of when the asymptotic solution develops 
remained open. The present author previously solved numerically similar equations for a 
quasi-one-dimensional conductor in a high magnetic field and found phase transitions to the 
DW state in both cases of repulsive and attractive interaction between electrons Q]. Even 
if initially the interaction was attractive (favorable for superconductivity), in the course of 
renormalization the sign of the interaction was effectively changed giving rise to a non-trivial 
DW state exhibiting the quantum Hall effect In present paper, parquet equations for an 
isotropic 3D metal in a very high magnetic field, with only the lowest spin-polarized Landau 
level filled, are solved numerically. It is found that in the repulsive case a phase transi- 
tion to the DW state occurs, in agreement with the analysis given in Ref. [Q]. In the case 
of attraction, no phase transition occurs and the system remains a non-trivially correlated 
metal to arbitrarily low temperatures. With decreasing temperature T, the system becomes 
progressively more one-dimensional, with the characteristic interaction length perpendicular 
to the magnetic field decreasing as C~^^^^ where ( = \n{uc/T), and is a cutoff of the order 
of the cyclotron frequency. However, the system never becomes strictly ID. By this reason, 
the SC susceptibility grows as (^^/'^ exp(constC^/^). This dependence is neither characteristic 
of ID, nor of conventional 3D metals. Thus, the system is a new non-trivial example of an 
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unconventional Fermi- liquid, neither of the Landau, nor of the Luttinger type . It should 
be emphasized also that there exists a model where the competition between SC and DW 
channels is suppressed and, thus superconductivity indeed may survive in an arbitrarily high 
magnetic field [|^. However, this is not the case for the model considered below. 

Let us consider an isotropic 3D metal with an effective mass m and a parabolic dispersion 
law. We assume that magnetic field is so high that only the lowest Landau level is filled 
and all electrons are spin-polarized. Dispersion law has the form e = p1/2m, where pz is 
the momentum along the field. The Fermi surface consists of the two Fermi points pz = 
±Pf = ±27f'^hnl'jj, where n is the volume concentration of electrons and Ih = (hc/eHY^'^ 
is the magnetic length. Coordinates X, Y in the plane perpendicular to the field will be 
measured in units of Ih- The wave functions of the lowest Landau level are 0: 

^.(X, Y) = n-'/%'/'exp{zxY - {x - X)V2). (1) 

Another representation can be obtained making a superposition of states (|l|): 



ijy{X,Y) = (27r)-i/2y exp{-txy)^,{X,Y) 

= 7,-'/%'" exp(zX(r -y)-{y- Yf/2). (2) 

To be distinguished from the running variables X and Y , the quantum numbers x and y, 
which label the wave functions (|l|) and (H), are denoted everywhere by the lowercase letters. 

Let us introduce now the operators a^{x,pz) and b'^{y,pz) which create electrons in the 
eigenstates (|1]) and (|^) with pz close to +pp and —pr, respectively. Near the Fermi energy, 
the electron spectrum can be linearized in pz and the Hamiltonian of the model can be 
written in the form H = Hq + Hi, where 

Ho = 2nlH J ^^^'P^^P^ ~ Pp") J dxa~^{x,pz)a{x,pz) 

-{Pz+Pf) J dyb^{y,Pz)b{y,Pz)], VF=PF/m, (3) 
^1= J 9-^-j00-^dx dy dx'dy' exp{-ixy' + iyx')'yo{x -x',y- y') 

d~^{x,pi)b+{y,P2)b{y',ps)d{x',pi+p2 -ps), (4) 



7o(r-r') = exp(zr Ar') J dX dY rx{X,Y)ry{X,Y)tljy,{X,Y)tlj,,{X,Y) 

= exp(-(r-r')V2). (5) 

In Eq. (1^), r = {x, y), r' = {x', y') and r A r' = xy' — yx' . The function 70 is a form-factor of 
the interaction between electrons in the representation given in (|1]) and (^. The interaction 
amphtude is (? = g2 — gi, where g2 and gi are the amphtudes of forward and backward 
scattering. We assume that the original interaction is not completely local, so that gi 7^ g2 
and g 0] however, its spatial range is much shorter than l^- 

The so-called parquet diagrams, which consist of the electron-hole and electron-electron 
loops inserted into each other in all possible ways, are the most important many-body 
corrections to the interaction vertex 70 (r) 0. The corrections form a series in powers of 

^ = {\g\/{27r)\FlH)H^c/max{T,uj,VF\\Pz\ - Pf\)], 

where a; is a frequency, since both of the one-loop diagrams are logarithmically divergent. 
The renormalized vertex of interaction 7(r, ^) obeys the following equation, shown graphi- 
cally in Fig.l: 

97(r,0/5e = / rfV7(r',07(r-r',0(l-e^^"'''), (6) 
7(r,0) = sign(^)7o(r). (7) 

Eq. (P) can be obtained from Eq. (2) of Ref. ^ via a Fourier transformation over the 
variable kx- 

The r.h.s. of Eq. is a difference of two terms. The first term is the contribution of the 
electron-hole loop, the second term — of the electron-electron loop. If only one of these terms 
is retained in the r.h.s., that corresponds to a ladder approximation instead of a parquet one, 
then Eq. (^) can be solved analytically. If one neglects the second term, then the equations 
are diagonalized via two-dimensional (2D) Fourier transformation 7(r,^) r(k, [Q: 

r(k,o = (sign(<7)ro-i(k)-o-\ (8) 
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where ro(k) is the 2D Fourier transform of 7o(r) (^. In the case g > (repulsion), solution 
(|) is called a "moving pole" 0, because the position ^p(k) = Fq ^(k) of the pole singularity 
in ^ depends on the value of k. It follows from that the minimum value of 'Cp(k), which is 
equal to = l/2vr, is attained at kc = 0. As was shown by a calculation of the appropriate 
susceptibility [Q, the moving pole singularity indicates a phase transition to a DW state, 
where the densities of the electron charge and spin are modulated along the magnetic field 
with a wave vector 2pF/h, and are homogeneous in the perpendicular plane (because kc = 0). 
Once the moving pole (||) develops, the second term in Eq. (||) indeed can be neglected, 
because this term contains an integration over k which makes it less singular than the first 
term 0. 

Solving numerically the full equation (P) aX g > 0, we find the moving pole singularity 
occuring at k^ = and = l-3/27r. The value of C,c is related to the transition temperature 
by the formula Tc = uJcexp{—{2n)^^cVFiH/\g\)- Thus, the only effect of the second term in 
Eq. (^ in the repulsive case is a certain decrease of the transition temperature. Otherwise, 
the ladder approximation gives qualitatively correct results. 

If the first term in the r.h.s. of Eq. (^) is neglected, then the resulting equation can be 
solved by performing the Fourier transformation 7(^1, y,^) A(xi,X2,0 o'^^^ ^^e variable 
y and introducing the function h{xi,X2,^) = A(xi — X2,xi + X2,0 obeying the following 
equation Q: 

dh{xi,X2,^)/d^ = — J dx h{xi,x,^)h{x,X2,^). (9) 
With initial conditions (0), Eq. (|^) has a solution: 

h{xi, X2, = V27cexp{-xl - X2)/(sign(^) + <). (10) 

In the case g < (attraction), Eq. (|T0|) has a pole singularity at = ^^at indicates 
a phase transition to a SC state. Solution (^) is called a "standing pole" 0, because the 
position of the pole in ^ does not depend on any continuous variable. For this reason, when 
expression (|To|) is substituted into the full Eq. (^, the first term has the same singularity as 
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the second, thus the SC ladder approximation can never be justified This fact explains 
why it is important to solve numerically the full Eq. (|^) in the case g < 0. 

Since the initial vertex 7o(r) (|^) depends only on r, which is the absolute value of r, then 
the same holds for 7(r,,^). Eq. can be rewritten for a new function 7(r, = 7(r, ,^), 
which depends on one spatial argument: 

roc rr+ri 



Jlr—ril 



sm 



[4(rir2)2-(r2 + ri-r2)2]i/2 ' 
7(r,0) = sign((7)exp(-rV2). (12) 



The numerical solution of Eq. (|TT])-(|T3) ^oi the case g < is shown in Fig. 2 for several 
values of "time" C,- After a short initial evolution, the function 7(r, attains the form 
7c(w(^)r), where u'(^) is a monotonically growing function of ^. The ansatz 7(r, ^) = 
7c(w(0^) is consistent with Eq. (|TT]) provided w is sufficiently large. In this case, the sine 
in Eq. ( p^D can be replaced by its argument, and Eq. (|1^) decouples into two equations: 

dw{0/d^ = Aw-\0, (13) 

foo rp+pi 
2Apd-^c{p)/dp = pidpi P2dp2lc{pi)lc{p2) 

■Jo Ap-Pi\ 

mpiP2r-{pi+pi-p')r', (14) 

A is a constant. It follows from (p!3[) that w{^) = [A{^ — ^o)]^''^- This dependence indeed 
was found numerically with = 0- The function 7c(p) is also known numerically: with the 
convention A = 1, 'jdp) = 7(p^~^^^5 0^ "where for the latter function one can take any of the 
plots in Fig. 2, except one corresponding to ,^ = 0. In summary, the solution of the parquet 
equations for the case g < has the self-similar form 

lir,0 = lcire^% (15) 

where we introduced a function of two variables 7c(r) = 'ydf). Eq. ( pTSD is neither a 
moving, nor a standing pole; it is rather a "squeezing" solution: the effective range of 
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interaction varies as ln~^^^(co'c/T), thereby making the system increasingly one-dimensional 
as temperature is reduced. 

In order to calculate susceptibilities, let us add to the Hamiltonian (|D-(|D the fictitious 
external fields, fsc and fow, which create electron-electron and electron- hole pairs: 



H2 = / —dxdy[fscix,y)a^{x,p;,)b^{y,-p^) + 
J Zn 

fDw{x,y) exp{-ixy)d~^ {x,p;,)b{y,p, - 2pf)] + h.c. (16) 

Let us start with the SC susceptibility. According to the parquet rules 0, it is necessary to 
calculate first a vertex |/,^), which is determined by the graphical equation shown in 
Fig. 3a: 

d^ir,0/d^ = - J rfr'7(r-r',0*(r',0e-^'-^^', (17) 
^(r,0) = /5c(r). (18) 

To solve this equation let us make a Fourier transformation over the variable y: ?/, ^) — >■ 

A{xi,X2,^) and fsci^i^v) ~^ Fsc{,Xi,X2)] then, introduce the new variables: x,,^) = 
A(a; — x, X + X, ^) and <l>5(7(x, x) = Fsc{x — x, x + x). In the new variables, Eq. ([T^) reads: 

5H(x, X, = - j dx' H(x, x', i)h{x\ X, 0, (19) 

S(x,x,0) = $5c(a:,x), (20) 

where the function h was introduced earlier. Eq. (p!9| ) is diagonal and degenerate with 

respect to x, the center-of-mass coordinate of a Cooper pair, so this variable can be omitted. 

Taking into account Eq. ([I5|) , one can rewrite Eq. ([T9|) in the form: 

c?S(x,^) f \ { ' ^' ^' c\ 

^ ~ J '^T/s'^'^y^ ' |T76 + '^T/3>^^^ + ^oTe ' 

where Ac(xi,X2) is the Fourier transform of 7c(xi,?/) over the variable y. An approximate 
solution is: 

^(x, = $5c(x) exp( — ), (21) 



Fo = - / 2Txr^c{r)dr = 3.1, 
Jo 

Fi = — / 7i:r^'yJr)dr = 1.1. 
Jo 



The contribution to the free energy 5iF due to the external field fsc is shown graphically 
in Fig. 3b: 

S,F = 2 f^dC [ dx\E{x,0\V\9\. (22) 







Substituting Eq. (^) into Eq. (p2D, one finds the susceptibility with respect to creation of 
an electron pair with the relative distance between electrons equal to x: 

X{x,0 = e^'exp{3Toe^' - T^x'e^y3)/\g\To. (23) 

When ^ — i> oo, the susceptibility diverges following an unusual law ^^/^ exp(const^^/'^). This 
behavior can be understood qualitatively in the following way. The total amplitude of 
scattering to all possible channels, 

7,(0 = Jdh^{v,o = -ro/e^', 

goes to zero with increasing ^. However, the function —2 J^'ytiOdC, which appears in the 
exponent of (|23|), increases as At the same time, the characteristic distance between 
electrons in a Cooper pair in (p3D squeezes as which indicates the increasingly one- 

dimensional character of the system. 

Following an analogous procedure, one can find a contribution to the free energy due to 
fnw- 



S2F= [^f^d^k\FDW 



1^1 

where FDw(k.) and Tc(k) are the 2D Fourier transforms of /dvk(i') and 7c(i")- For all values 
of k the DW susceptibilities have finite limits at zero temperature. 

The lowest order diagram, which gives a logarithmic correction to the one-particle Green 
function, is shown in Fig. 4. Renormalized Green functions G±, where the indices ± refer 
to the d and b fermions respectively, can be found in the framework of a scaling hypothesis 
13: 



G,fc..) = ^^Et2£HqW4(2^_ (24) 

ujtvf{PzTPf) 



roc 



Fa = / 2TTr%{r)dr = 1.9. 
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Eq. (p^) differs from expressions known for tlie 3D Landau and tlie ID Luttinger Fermi- 
liquids 0. 

In conclusion, the self-similar "squeezing" solution ( pTSD of the parquet equations has 
never appeared previously in condensed matter physics. The correlation functions (|23| ) and 
(P^) represent a new universality class of Fermi-liquid behavior. All of this enriches our 
intuition beyond the Luttinger and the Landau schemes in consideration of what happens 
when many ID channels interact []Tl |. 

The author is grateful to E. Abrahams for the support of this work via the NSF Grant 
No. DMR 89-06958. 
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FIGURES 

FIG. 1. Parquet equations for the vertex of interaction 7(^2 — xi, y2 — yi, shown as a filled 
circle. The bare vertex ^0(^2 — xi,y2 — yi) is denoted as a dot. The Green functions of electrons 
with momenta pz close to ±pF are shown as solid and dotted lines. 

FIG. 2. The vertex of interaction 7(r, ^) as a function of r at 27r^ = (curve a), 5 (b), 500 (c) 
and 3200 (d). 

FIG. 3. a) Parquet equations for the vertex shown as a filled triangle. The external 

field fscix, y) is denoted as a square, b) Contribution to the free energy 5iF due to the field fsc- 



FIG. 4. The lowest order logarithmic correction to a Green function. 
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